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Abstract 

We investigate the formation of light nuclei with the nuclear mass number less than or equal 
to four in 2+1 flavor QCD using a non-perturbative improved Wilson quark and Iwasaki gauge 
actions. The quark mass is decreased from our previous work to the one corresponding to the pion 
mass of 0.30 GeV. In each multi-nucleon channel, the energy shift of the ground state relative to 
the assembly of free nucleons is calculated on two volumes, whose spatial extents are 4.3 fm and 
5.8 fm. From the volume dependence of the energy shift, we distinguish a bound state of multi 
nucleons from an attractive scattering state. We find that all the ground states measured in this 
calculation are bound states. As in the previous studies at larger m n , our result indicates that at 
m.„- = 0.30 GeV the effective interaction between nucleons in the light nuclei is relatively stronger 
than the one in nature, since the results for the binding energies are larger than the experimental 
values and a bound state appears in the dineutron channel, which is not observed in experiment. 
Possible sources of systematic error in this calculation are discussed. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38.-t 12.38.Gc 
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I. INTRODUCTION 


The strong interaction is the origin of the formation of nuclei. Non-perturbative lattice 
QCD calculation is a powerful tool to confirm nucleus formation from the first principle of 
the strong interaction. The nucleus formation was examined in lattice QCD in Ref. 


A 


m 


which the binding energies for the 4 He and 3 He nuclei were calculated in quenched QCD at 
a heavy quark mass corresponding to the pion mass m n = 0.80 GeV. In this calculation, 
a multi-nucleon bound state was identified by the volume dependence of energy shift of 
the ground state relative to the assembly of free nucleons. This study was followed by 


calculations in Nf = 3 QCD at m n = 0.81 GeV [2] and 2 + 1 QCD |3| at m n = 0.51 GeV 
The 4 He nucleus formation was also reported in a different approach usingjdie two-nucleon 
potential calculation in Nf = 3 QCD at m n = 0.47 


energy reported is much smaller than those in Refs, jl 


.02, and 1.17 GeV |4|. The binding 


3]. The authors in Ref. [1] also made 


the first systematic study of the volume dependence of the energy shift for the two-nucleon 


states in the spin triplet 3 Si and singlet 1 Sn channels in quenched QCD [5]. This work was 


extended to the full QCD case in Refs. 


the ear 


ier studies of these channels 


0 


? Is, 6]; The volume dependence was not examined in 
|8|. The calculation in Nf = 2 + 1 QCD at m n = 0.39 
GeV |6| was not conclusive of bound state formations in both channels due to large errors 

n n rt 

of the energy shift. Other calculations [2|, |3|, 5], on the other hand, concluded that there 

lannel. The latter results conflict with the one from the two-nucleon 

9 |- 


is a nucleus in each c 
potential calculation 


The results obtained from calculation of the energy shift seem to indicate that the effective 
interaction among nucleons seems relatively stronger, compared to the kinetic energy of the 
nucleons, than those in nature. Indeed, in the calculations done to date, the binding energies 
for the 3 He and deuteron are clearly larger than the experimental values; also there is a 
bound state in the dineutron channel ( 4 So channel), which has not been observed in nature. 
A possible explanation of the discrepancy between the lattice QCD results and experiment is 
the heavier u, d quark masses employed in the calculations. If this is the case, at the physical 
quark mass the binding energies would agree with those in the nature, and the bound state 
in the 4 So channel would disappear. In order to check this scenario, calculations at lighter 
quark masses than those employed in the previous calculations are necessary. In this work, 
we extend our previous calculation at m n = 0.51 GeV in 2+1 flavor QCD to a smaller 
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quark mass of m n = 0.30 GeV. We investigate whether the light nuclei ( 4 He, 3 He, deuteron, 
and dineutron), which were observed in the previous calculations [l-l], 5], are formed or not 
at this quark mass. 

This paper is organized as follows. In Sec. [TT] we explain details of calculation including 
the simulation parameters for gauge configuration generation and the calculation method 
for the multi-nucleon channels. Section m presents the results for the 4 He, 3 He, deuteron 
( 3 Si), and dineutron ( 1 S 0 ) channels. Comparisons of our results with those in the previous 
studies are also discussed. Conclusions and discussions are given in Sec. HV1 


II. SIMULATION DETAILS 


A. Simulation parameters 


For 


action 


gauge configuration generation in 2+1 flavor QCD, we employ the Iwasaki gauge 


10|] and a non-perturbative 0(a)-improved Wilson quark action. The bare coupling 
is fixed ai (3 = 1.90 for which we use csw = 1-715 fill ]. The lattice spacing is a = 0.08995(40) 
fm, corresponding to a -1 = 2.194(10) GeV, which was determined by tuq = 1.6725 GeV 12]. 


We choose two lattice sizes, L 3 x T = 48 3 x 48 and 64 3 x 64, to investigate the spatial 
volume dependence of the energy shift between the multi-nucleon ground state and the free 
nucleons. The physical spatial extents are 4.3 and 5.8 fm, respectively. We choose the 
hopping parameters («w, ac s ) = (0.1376314,0.1367299) to obtain m n = 0.30 GeV and the 



marize the simulation parameters in Table H including the block sizes and the preconditioning 
factor in MPDDHMC and the polynomial order in UVPHMC. We take r = 1 for the tra¬ 
jectory length of the molecular dynamics in all the runs. The step sizes are chosen such 
that we obtain the reasonable acceptance rates presented in Table [T] In the spatial extents 
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of 4.3 fm and 5.8 fm, 1000 and 800 trajectories are generated in four and two streams after 
thermalization, and the total lengths of trajectory for the measurement are 4000 and 1600, 
respectively. 

We calculate correlation functions in the multi-nucleon channels in every 10 trajectories 
for both volumes using the same quark action as for the configuration generation. The er¬ 
rors are estimated by jackknife analysis choosing 200 and 160 trajectories for the bin size for 
the smaller and larger volumes, respectively. Statistics is increased by repeating the mea¬ 
surement of the correlation functions with different source positions on each configuration. 
We calculate the correlation functions not only in the temporal direction but also in the 
spatial ones with the use of the space-time rotational symmetry. It allows to increase the 
statistics by a factor four effectively. The parameters of the measurement, e.g., the number 
of configurations and the bin sizes, are listed in Table [Til 

B. Calculation method 

We extract the ground state energy in the multi-nucleon channels and the nucleon mass 
from the correlation function, 

Go(t) — (0|(9(f)(9(0)|0), (1) 

with O being proper operators for the 4 He, 3 He, 3 Si and 1 So channels and also the nucleon 
N, which are given in the next subsection. 

We define the energy shift between the multi-nucleon ground state and free nucleons on 
finite volume as, 

AEl = Eq — N^rnjy, ( 2 ) 

where Eq is the lowest energy level for a multi-nucleon channel, Nn is the number of nucleons 
in the channel, and is the nucleon mass. This quantity is directly extracted from the 
ratio of the multi-nucleon correlation function to the iVvth power of the nucleon correlation 
function 

Gp{t) 

(< G N (t )) N "’ 

in the large time region where both correlation functions are dominated by the ground state. 
We also define an effective energy shift as 

A Ef = In ( ) , (4) 
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which is utilized to investigate plateau region in later section. Note that the definition of 
A Ei and A E e ^ follows that in Refs, js, 5], while the sign convention is opposite to that in 
Ref. fl. 

There are two computational difficulties in the calculation of Go(t) for multi-nucleon 
channels. One is a factorially large number of Wick contractions for quark-antiquark fields. 
To overcome the difficulty, we use the reduction technique of calculation cost proposed in 
our exploratory work [l|. It is noted that other reduction tec’ 
of Wick contractions have been proposed for the multi-meson 


3 and multi-baryon 


num 

ber 

21- 

-23] 


channels. Another problem is an exponential increase of statistical errors with atomic mass 
number. For this difficulty, we carry out measurements as much as possible using multiple 


actor twelve and five times larger than 
for 4.3 fin and 5.8 fm spatial extents, 


source points. The number of measurements are a 
those in the previous calculation of = 0.51 GeV 
respectively. 

Another difficulty in the nucleus calculation is to distin guis h a multi-nucleon bound state 
from an attractive scattering state in a finite volume 24d26l]. This problem is handled by 
studying the volume dependence of the measured A El as in Refs, [l|, [5]. While the en ergy 
shift of an attractive scattering state vanishes in the infinite volume limit as 1/L 3 24 1271]. 
the physical binding energy of a bound state remains at a finite value in the limit. 


C. Interpolating operators 

The u, d quark propagators are solved with the periodic boundary condition in all of 
spatial and temporal directions using an exponentially smeared source, 

V(x,t) = ^2Ae- BlS ^q 0 (y,t), (5) 

y 

for |£| 4 0; and q(x,t) = 1 for |x| = 0, after the Coulomb gauge fixing, where q 0 is 
the local quark field. We choose the smearing parameters A = 0.8 and B = 0.16 on 
the two volumes to obtain reasonable plateaus of the effective energy for the nucleon and 
ground states in the multi-nucleon channels. The stopping condition of the quark propagator 
e = |Dx — b\/\b\ < 10” 6 is applied in both volumes to reduce the calculation time. We have 
checked in a subset of the configurations that the results with this looser stopping condition 
agree with the ones using a more stringent stopping condition e < 10” 11 in more than six 
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digits, even for Gh He (t) at t — 12. The systematic error coming from the discrepancy is 
much smaller than the statistical error in the current calculation. 

The interpolating operator for the proton is defined as 

Pa = £abc([Ua] t C'y 5 d b )u “, ( 6 ) 


where C = 7472 and a and a, b, c are the Dirac and color indices, respectively. The neutron 
operator n a is obtained by replacing by d" in Eq. (J 6 ]). 

The 4 He nucleus has zero total angular momentum and positive parity J p = 0 + with the 
isospin 7 = 0. We employ the simplest 4 He interpolating operator with zero orbital angular 


momentum L — 0, and hence J = S with S being the total spin 


28|, 


where 


4 He 


1 

72 


(x:n - xv), 


( 7 ) 


x - — 1 —] + H - H-H _ [ _ + -I—D> 


( 8 ) 


X ~ , ^=([d-1—] + [- H-h] + [H-f] + [- + H—] - 2[+ H-] - 2[-1—h]),( 9 ) 

with +/— being up/down spin of each nucleon, and 77,77 are obtained by replacing +/— in 
X,X by p/n for the isospin. 

The 3 He nucleus has J p = ^ + , I = \ and I z = We employ the interpolating operator 


in Ref. 


m, 


3 He = (\p-n + p + ) - \p + n+p_) + \n + p + p_) - \n + p_p + ) + \p + p_n + ) - \p_p + n + )). ( 10 ) 


The two-nucleon operators for the 3 Si and 4 So channels are given by 


NN 3Sl (t) = 7= [p + (t)n + (t] - n+(t)p+(t )], (11) 

NNi So (t) = 7= \p+(t)p_(t) - p-{t)p + {t )\, (12) 

respectively. In the 3 Si channel the operators for the other two spin components are con¬ 
structed in a similar way. We increase statistics by averaging over the three correlation 
function with each spin component operator. 

Using the interpolating operators above, we calculate correlation functions in each chan¬ 
nel. For the source operators in all correlation functions, we insert the smeared quark fields 
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of Eq. (J5J) for each nucleon operator located at the same spatial point x. Each nucleon in 
the sink operator, on the other hand, is composed of the point quark fields corresponding 
to qo in Eq. (J5]), and projected to zero spatial momentum. To save the computational cost 
we use the non-relativistic quark operators, in which the Dirac index in Eq. fl6]) is restricted 
to the upper two components in the Dirac representation. 

III. RESULTS 

A. Nucleon and pion masses 

The results for effective m n in the two volumes are shown in Fig. [Tj together with the 
exponential fit result of Cjv(t) and the one standard deviation error band. Plateaus are 
clearly seen for t > 8 for both volumes. The difference of the fit results between the two 
volumes is 1.4 standard deviations, and hence statistically not very significant. We also 
do not expect much finite size effect for these large volumes satisfying m n L >6. In the 
following sections, we therefore consider that the difference is caused by statistics, and will 
not estimate the systematic error from it. The pion effective masses in each volume show 
better consistency than the nucleon mass, as presented in Fig. [2j Those fit results are 
tabulated in Table El 

B. 4 He channel 

The effective energy energy shift Adefined in Eq. 03]) is shown in Fig. [H] for the two 
volumes. Clear signals are seen for t < 10, but for larger t the statistical error increases 
rapidly. A plateau appears at t = 9-12 on the larger volume, while it is not clearly seen 
on the smaller volume. More statistics is desirable for establishing a plateau in this case. 
We calculate the energy shift A El in Eq. 02]) by a single exponential fit of R(t) in Eq. 02]) 
using the same range t = 9-12 for the two volumes. The systematic error is estimated from 
the variation of the fit results with six different fit ranges, where the minimum or maximum 
time slice is changed by ±1, and the minimum and maximum time slices are changed by +1 
and +2. The central fit result is shown in Fig. [3] by solid lines with the band representing 
the statistical error. The dashed lines denote the total error adding the statistical and 
systematic errors by quadrature. In Fig. [4] we illustrate how we estimate the systematic 
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error. Shown in the figure are the results of 8 fits obtained by shifting the fitting range as 
explained in the figure caption. The horizontal band with solid lines shows the total error 
obtained by adding the statistical error and the systematic error from the 7 fits on the left 
by quadrature. We observe that the 7 fits on the left reasonably covers the variation, with 
the 8th fit at the rightmost with the fitting range shifted by +3 from (t min ,t max ) falling 
within the band of solid lines within one sigma. We therefore consider that our estimate 
of systematic errors is reasonable under the current statistics. The values of A El with the 
statistical and systematic errors are summarized in Table IIII1 

Figure [5] shows the volume dependence of A El as a function of 1/L 3 . The inner bar 
of each data denotes the statistical error and the outer bar represents the total error with 
the statistical and systematic ones added in quadrature. Since the volume dependence is 
not large, we estimate the energy shift in the infinite volume limit AE^ by a constant fit 
as presented by solid line and open circle in Fig. [5j An exponential type extrapolation, 
A El = AEqq + C exp (—C e L), cannot be carried out in this work due to smaller number of 
data than its free parameters. The systematic error is estimated from the variation of the 
central values obtained by 49 fits. The 49 fits are constant fits with various combination of 
7 2 data set, where in each volume we have 7 data with different fit range of R(t ) as explained 
in the above. The result of A FA with the statistical and systematic errors are tabulated in 
Table eld From the result that A FA, is non-zero and negative, we conclude that the ground 
state is bound in this channel. The binding energy equals —A = 47(7) MeV where 
the first and second errors are statistical and systematic, respectively. 


The result for —AAV is compared with the experimental value of 28.3 MeV and with the 
previous three results [1-3] in Fig. [HI The binding energy for m n = 0.3 GeV obtained in this 
work is similar in magnitude with our previous results for Nf = 2 + 1 at m n = 0.51 GeV 


and Nf = 0 at m n = 0.80 GeV [5]. Compared to experiment, if one used the upper total 
error, our current value is consistent within 1.5 a. The result of the Nf = 3 calculation at 


m 7 r = 0.81 GeV [2] is about three times larger than the other results. This difference might 
be due to different quark masses of the calculation or dynamical quark effects. On the other 
hand, the result obtained with the two-nucleon potential extracted from Nf = 3 calculations 
at m n = 0.47 GeV [4] has a very small binding energy, A E = 5.1 MeV, compared to the 
other results. 






C. 3 He channel 


Figure [3 shows the effective energy shift Ain Eq. (J4]) for the two volumes. The signals 
are better than those in the 4 He channel shown in Fig. [3l A plateau is seen for the smaller 
volume case, while it is less clear in the region of t — 8-12 for the larger volume case. The 
energy shift A El in Eq. (J2J is determined by an exponential fit to R(t) in Eq. (J3]) with 
the fit range of t — 8-12 and 9-12 for the smaller and larger volumes, respectively. The 
systematic error of A El is estimated in the same way as for the 4 He case as described in 
the above subsection. The fit results with the statistical and systematic errors are shown 
in Fig. [T| and Table IIIII The explanations for the solid and dashed lines are given in the 
previous subsection. Figure [8] shows how we estimate the systematic error. The relative 
difference of the extra fit result with the slid fit range by +3, as explained in the previous 
subsection, from the central fit result with the total error is less than 1.6 a. 


A wea 
results 


;ci 

0 


c volume dependence of A El observed in Fig. [9] is similar to those in the previous 


3]. A constant fit of A El with only the statistical error gives a large value of 
X 2 /d.o.f. = 4.1. It agrees with the two data within the total error as shown in Fig. EE 
however. Thus, we take the constant fit result as the estimate of the central value of AE^ 
in this calculation. 

We estimate the systematic error of A E^ in the same way to the 4 He case. We omit 
18 constant fit results, however, with y 2 /d.o.f. > 4.1. The extrapolated result of Ais 
clearly nonzero and negative as presented in Fig. EH Thus the ground state is a bound state, 
corresponding to the 3 He nucleus, with a binding energy of —AE^ = 21.7(1.2)(t} 3 6 ) MeV, 
where the first and second errors are statistical and systematic, respectively. 

The quark mass dependence of the energy shift is plotted in Fig. [10J Our present re- 
suit together with our two prev.ous calculations Q S show very small dependence, while 
NPLQCD reported a much deeper bound state [2|. All lattice results in the figure have the 
binding energy larger than the experimental value 7.72 MeV. 


D. Two-nucleon channels 

We present AEy^ in Eq. (J4]) for the 3 Si and 4 So channels in Figs. []J]and[T2l respectively. 
The signals are clean up to t ~ 14, but statistical fluctuations spoil the signals in the 
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in the 1 So channe 


ire 

previous studies [2 

•• 

m 


are_smaller than those in the 3 Si 
We observe a clear plateau 
with a negative energy shift for 9<7<^14, although of a less quality for the 3 Si channel for 
the smaller volume. We determine A El by an exponential fit to R(t) of Eq. OSJ) with the 
fixed fit range of t — 9-13 for the 3 Si channel, and with t — 10 — 14 for 1 So- The fit results 
are presented in Figs. |TT] and |T21 and are summarized in Table IIVI The systematic error 
estimations for the 3 Si and 1 So channels using the results with the several fit ranges, as in 
the 4 He and 3 He cases, are presented in Figs. H3]and[l4l respectively. 

The volume dependence of the energy shift in the two channels is shown in Figs. flSlandflbl 
In both channels, the volume dependences arc mild, so that the data can be reasonably fitted 
by a constant. The results for the constant fit are non-zero and negative. This indicates 
that the ground states in the two channels are bound states. The same conclusion is also 


obtained by a fit including finite volume effects of the two-particle bound state 

.,2 


25 


26|, 


A E t . = 


7 


m N 


1 + ^2 ^ exp(-yLVrT) \ 


1 L 




(13) 


where 7 and C 7 are free parameters, n is a three-dimensional integer vector, and denotes 
the summation without \n\ = 0. In the fit, we use the weighted average value of mjv with 
the two volume data. In the above equation, it is assumed that 


A Eqo — 


7 

m N 


2 m 


N 


m 


N 


7 


(14) 


Note that the degrees of freedom are zero in the fit with Eq. (fT3]h The fit result is presented 
in each figure at 1/L 3 = 0. We take the constant fit as the central value of the binding 
energy —AEoo, and estimate the systematic error in the same way as in other channels. In 
the systematic error estimation, we include the fit result using Eq. (113]) . while we exclude 9 
and 8 constant fit results in the 3 Si and 1 So channels, respectively, which yield y 2 /d.o.f. > 3. 
The results for the binding energy are —AEoo = 14.5(0.7) (i^g) MeV for the 3 Si channel 
and 8.5(0.7)(1JJ) MeV for 1 S 0 with the first and second errors being the statistical and 
systematic, which are also summarized in Table IIVI 

In Figs. H7]and[IEl the results for AE^ in the present work are compared with those of 
the previous studies [2, 3 i5HS] as a function of m 2 . Almost all results report negative values, 


except for those of Ref. [8] with large errors. The earlier calculations |7|,|8| did not investigate 


the volume dependence of A El. More recent studies 


[2, a, i Q 


fll 


examined the dependence 
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and estimated the infinite volume value through extrapolations 


there is no significant volume dependence of A Ei 


0 - 


b y, 3 


or checked that 


All the recent results suggest that 
the ground states in both the channels are bound states. One exception is Ref. j(3] where the 
conclusion is not clear due to large errors. 

While lattice results are mutually qualitatively consistent, they differ from experiment 
in more than one aspects. For the 3 Si channel, the binding energy —Afound in the 
lattice calculations [2, 3, 5, 6] is a factor five to ten times larger than the experimental value. 
Furthermore, we observe no tendency in the binding energy to approach the experimental 
value, at least over the pion mass range ro, = 0.3-0.51 GeV. For the 1 So channel, the bound 
state found in the lattice calculations is absent in experiment. Furthermore, similarly to the 
3 Si channel, the binding energy is almost flat in m\ in the interval m n = 0.30-0.51 GeV. It 
is not clear whether the bound state observed in the lattice calculation becomes unbound 
toward the physical m 7T . 


IV. CONCLUSION AND DISCUSSION 


We have extended our previous nuclei calculation in 2+1 flavor QCD at = 0.51 
GeV j^] to the lighter quark mass corresponding to m n = 0.30 GeV and m N = 1.05 GeV. 
In order to suppress an exponential increase of statistical errors at smaller m„ we have 
carried out a much larger number of measurements by a factor twelve and five for the case 
of the spatial extent of 4.3 fin (48 3 ) and 5.8 fin (64 3 ), respectively, compared to those for 
the m n = 0.51 GeV case with the same volumes. We have found that in all channels we 
have studied, 4 He, 3 He, and two-nucleon 3 Si and 1 So, the ground state is a bound state by 
investigating the volume dependence of energy shift A El. The binding energies estimated 
for the infinite volume are as follows: 

47(7) (if?) MeV for 4 He, 

21.7(1.2) (+} 3 6 ) MeV for 3 He, 

14.5(0.7)(l^; 4 ) MeV for 3 Sx, 

8.5(0.7)MeV for 4 S 0 . 

These values differ little from those obtained at m n = 0.51 GeV j^|. The largest relative 
difference occurs for the 3 Si channel, which is only a 1.9 a effect if we use the total error 
adding the statistical and systematic ones by quadrature. Therefore, our conclusions at 


- AFW = 


(15) 
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= 0.30 GeV are similar to those in Ref. [3| for m n = 0.51 GeV: the binding energy of 
the 4 He nucleus is comparable with the experimental value, while the 3 He nucleus and the 
deuteron are about three and seven times larger than the experimental values, respectively, 
and a bound dineutron is observed in the 4 So channel. 


The differences we observe from experiment may arise from various sources, either com¬ 
putational or physical in origin. Statistical errors are fairly large in the calculations even 
for light nuclei. While the negative value of the energy shift is certain in all channels we 
looked at, better statistics and improved techniques will be welcome to better control the 
extraction of the energy shift for each volume and the infinite volume extrapolations. 

The quark mass is heavier than experiment in all calculations to date. The binding in the 
1 So channel is shallower than the 3 Si so that the former bound state may become unbound 
as decreases toward the physical value. This can only be verified by calculations of the 
nuclear binding energy at smaller quark masses. 


It is also possible that finite lattice spacing effect is rather subtle. The short distance 
repulsion, in the language of nuclear potential, is possibly affected more by such effects 
than the long distance attraction, so that finite lattice spacing effects may push out multi¬ 
nucleon wave function, and then the ground state would become a scattering state, for lattice 
spacings smaller than some value. 


Another possible source of systematic error is excited state contaminations in the calcu¬ 
lation. We have assumed that the nucleon and nucleus correlation functions are dominated 
by the ground state in the large t region, where the plateau of AE in Eq. (j4j) appears. 
While we have tuned the smearing parameter of the quark field to increase the overlap of the 
nucleus operator to its ground state, from the current data we cannot completely exclude 
the possibility that it is not sufficient to suppress the contaminations. To investigate the size 


of possible contaminations, we might try analyses with the variational method 
correlation function matrices. 


30| using 


For now, however, we think that a calculation at the physical point, keeping the lattice 
spacing, is the next step. 
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TABLE I. Simulation parameters for gauge configuration generation at («w, 
k s )=( 0.1376314,0.1367299). The definition of parameters is same as in Ref. [13]. 


L 3 x T 

48 3 x 48 

64 3 x 64 

# run 

4 

2 

(-No, N\ , IV 2 , IV 3 ) 

( 2 , 2 , 2 , 6 ) 

( 2 , 2 , 2 , 8 ) 

Block size 

12 2 x 6 2 

OO 

W 

X 

P 

0.998 

0.998 

Apoly 

320 

340 

MD time 

1000 1000 1000 1000 

800 800 

P acc (HMC) 

0.818 0.801 0.814 0.794 

0.902 0.880 

P acc ( GMP) 

0.959 0.959 0.962 0.962 

0.954 0.967 


TABLE II. Number of configurations, separation of trajectories between each measurement in the 
units of r, bin size in jackknife analysis in the units of configuration, number of measurements 
on each configuration, exponential smearing parameter set ( A,B ) in Eq. ([5]), pion mass m, and 
nucleon mass mjsr are summarized for each lattice size. Number of measurements include factor 
four by measurement with all four directions, which is explained in the text. 


L 

T 

# conhg. 

7"sep 

bin size # rneas. 

(A,B) 

m n [GeV] 

rriN [GeV] 

48 

48 

400 

10 

20 1152 

(0.8,0.16) 

0.3001(14) 

1.057(2) 

64 

64 

160 

10 

16 2048 

(0.8,0.16) 

0.2987(9) 

1.053(2) 


TABLE III. Energy shift —A El in physical units and fit range for ^He and 3 He channels on each 
spatial volume. Extrapolated results in the infinite spatial volume limit are also presented. The 
first and second errors are statistical and systematic, respectively. 



4 He 


3 He 


L 

-A E l [MeV] 

fit range 

-A E l [MeV] 

fit range 

48 

46(13)(±“) 

9-12 

18.7(1.9X1^) 

8-12 

64 

47(9)C|) 

9-12 

23.7(1.6)(1 43 5 ) 

9-12 

OO 

47(7)(1 2 ?) 

— 

21.7(1.2)(1} 3 ) 

— 
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TABLE IV. Same 

as Table Hill for 3 Si 

and 1 So channels. 



3 Si 


'So 


L 

-A E l [MeV] 

fit range 

-A E l [MeV] 

fit range 

48 

13.8(0.9)(±f;f) 

9-13 

7.7(0.9)(±?;|) 

9-13 

64 

15.6(1.2)(t};°) 

10-14 

9.5(0.9)(±g;i) 

10-14 

oo 

14.5(0.7)(±g;|) 

— 

8.5(0.7)(±S;g) 

— 
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FIG. 1. Nucleon effective masses on (4.3 fm) 3 and (5.8 fm) 3 volumes in lattice unites. Fit result 
with one standard deviation error band is expressed by solid lines. 



FIG. 2. Same figure as Fig. [U but for pion effective masses. 
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FIG. 3. Effective energy shift Afor 4 He channel on (4.3 fm) 3 (top) and (5.8 fm) 3 (bottom) 
volumes in lattice units. Fit result with one standard deviation error band and total error including 
the systematic one is expressed by solid and dashed lines, respectively. 
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FIG. 4. Fit range dependence of energy shift A El for 4 He channel on (4.3 fm) 3 (top) and 
(5.8 fm) 3 (bottom) volumes in lattice units. The horizontal axis corresponds to the fit range 

(^miru^max); (^min Ij^max); (^min + 1 5 ^max)i (^min? ^max l)j (^min? ^max +i), (t min +i 5 ^max +i), (t min + 

2,t max + 2), and (t m i n + 3,f max + 3) from left to right. f m i n and f max are minimum and maximum 
time slices of the fit range, respectively, whose values are explained in the text. The total error 
band including the statistical and systematic is expressed by solid lines. 
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FIG. 5. Spatial volume dependence of A El for 4 He channel in MeV units. Outer bar denotes the 
total error of statistical and systematic ones added in quadrature. Inner bar is for the statistical 
error. Constant fit result is shown by open circle symbol. Experimental value (star) is also 
presented. 
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FIG. 6. m ^ dependence of energy shift for ^He channel in MeV units. Open and closed symbols 


0 


BO 


0 ] 


denote the quenched [if] and full QCD [2|i4[] results, respectively. The results of Refs. Jl|, f3| and 
this work are the ones in the infinite volume limit. The error of the result obtained from the 
two-nucleon potential was not estimated in Ref. [J]. Experimental result (star) is also presented 


for comparison. 
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FIG. 7. Same figure as Fig. [3j but for 3 He channel. 
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FIG. 8. Same figure as Fig. [4J but for 3 He channel. 
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FIG. 9. Same figure as Fig. 0 but for ‘^He channel. 
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FIG. 10. Same figure as Fig. (6[ but for 3 He channel. 
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FIG. 11. Same figure as Fig. O but for 3 Si NN channel. 
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FIG. 12. Same figure as Fig. O but for 1 Sq NN channel. 
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FIG. 13. Same figure as Fig. SI but for 3 Si NN channel. 
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FIG. 14. Same figure as Fig. SI but for 1 Sq NN channel. 
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FIG. 15. Same figure as Fig.[5l but for 3 Si NN channel. Fit result using finite volume dependence 
of two-particle binding energy Ea. (fT3l) is also plotted. 
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FIG. 16. Same figure as Fig. [15l but for 1 Sq NN channel. 
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FIG 


. 17. Same figure as Fig. [6l but for 3 Si 


NN channel. Open and closed syrnbo 


s denote the 


quenched [5j, |7(] and full QCD |2|, yj, [y, [8( results, respectively. The results of Refs. [2 
this work are the ones in the infinite volume limit. 
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FIG. 18. Same figure as Fig. [T7I but for 1 Sq NN channel. 
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